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The overlap of a d + 1 dimensional directed polymer of length t in a random 
medium is studied using a Renormalization Group approach. In d > 2 it vanishes 
at Tc for t — > oo as t^ where S = [jEh\i and z is the transverse spatial rescaling 
exponent. The same formula holds in d = 1 for any finite temperature and it agrees 
with previous numerical simulations aX d = 1. Among other results we also determine 
the scaling exponent for mutual repulsion of two chains in the random medium. 



64.60.Ak, 05.40+J, 75.10 Nr, 36.20-r 



Typeset Using REVTEX 



1 



One commonly studied quantity in random systems is the overlap of appropriate physical 
variables in different states as, for example, the overlap of magnetizations in two states for a 
spin glass system |1|J^. The distribution function for the overlap q, p{q), which, e.g, in a pure 
Ising system has two 6 function peaks for the two possible low temperature ferromagnetic 
states, has a nontrivial structure in the spin glass phase because of the broken ergodicity. 
Such nontriviality of the overlap distribution is the main characteristic for spin glass type 
systems. Notwithstanding the importance of overlaps, very little is understood analytically 
about it. 

Of late, a directed polymer in a random environment is taken to be the paradigm of 
disordered systems. This is because of the strong analogy in behavior with other random 
systems and spin glass P], and more so because of the availability of exact results through 
Bethe ansatz Q], use of various analytical and numerical methods like, nonlinear differntial 
equation transfer matrix and other approaches. Even in this situation, the question 
of overlap remains a mind boggling issue 0J§]- In this paper, we implement a renormalization 
group (RG) approach to obtain the scaling behavior of the overlap. To our knowledge, no 
such analytical result for overlap is known for any other realistic random system. 

The d + 1 dimensional directed polymers are random walks stretched along a particular 
direction with fluctuations in the transverse d dimensional space. These string like objects 
appear in various random systems of interface fluctuations and pinning crystal growth 
1^, spin glasses etc. Most of the recent work done on this problem attempt to understand 
the low temperature (strong disorder) phase which, as a matter of fact, is the only possible 
phase in 1 + 1 dimensions. The problem in 1 + 1 dimensions is almost settled with relevant 
exponents known exactly |T0|. The two important exponents are x ^^^^ C = 1/^5 which 
describe the free energy fluctuation and transverse size as the length t ^ oo, namely / ~ t^^^ 
and < X >~ t''. For (i = 1, x = 1/2 and z = 3/2, with x + z = 2. However, precise values of 
these exponents in higher dimensions and their exact dimensional dependences are not yet 
well understood mainly because of the lack of any perturbative fixed point |TT]]. For high 



enough dimensions {d > 2) it is found that there is a phase transition [|l3,|T^ from the high 
temperature (weak disorder) to the low temperature phase (strong disorder). What happens 
at (i = 2 is not clear though the consensus seems to be against a phase transition 
We mention, in passing, that several exact results are known, especially in connection with 
such disorder induced phase transitions, if the randomness is in the interaction instead of 
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The high temperature phase of a directed polymer in a random potential in c? > 2 
is simple since there the quenched and annealed free energies are equal. The situation 
is more complex in the low temperature phase because of the nonzero overlap and the 
subsequent nonanalyticity of the free energy that supports the hypothesis of the coexistance 
of several pure states - a phenomena reminiscent of spin glasses. Partial information 
regarding the thermal, geometrical properties and a few related distribution functions in 
the low temperature phase are available from Monte Carlo simulations, expansion methods 
and transfer matrix techniques P^|8|. However, at present the overlap, especially near the 
critical point, seems to elude these techniques. In this paper our main focus is at this 
transition temperature. Unlike the tree problem which, in some way, corresponds to the 
mean field limit d — > oo, our result is true for finite dimensions. [See, e.g., Ref. [|16| for the 
peculiarities of the Cayley tree problem.] 
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Similar to the concept of the overlap of two different magnetization states here in the 
DP picture, the overlap means the number of common spatial points visited by two different 
configurations of the polymer. Introducing two different configurations of a polymer is 
equivalent to starting with the original polymer with a replica. The overlap is then the 
average number of contacts of these two polymers (see below for a more precise definition). 
The procedure for the evaluation of the overlap would be to introduce a new interaction 
that penalizes such contacts with coupling constant vq. The overlap follows, as statistical 
mechanics prescribes, from the calculation of the appropriate derivative of the free energy 
of such an interacting system. 

In the path integral representation the working hamiltonian for two interacting directed 
polymers 



t 



^ = /^^ [ E (ixi^W + ^nxi(r),r)j +-.;o5(xi2(r))] (1) 

where Xj(t) is the d dimensional spatial coordinate of the zth chain at the contour length 
t, Xj(t) = and Xi2 is the relative seperation of the two chains. The first two terms 

represent the entropic fluctuations of two free Gaussian chains with 7 as the bare line tension. 
V corresponds to the space and time dependent random potential seen by the two chains 
interacting at same t with 5 function potential of strength \vq/{2'^). The significance of 
the peculiar factor ^ with the random potential and with the coupling constant Vq will be 
clear from discussions later. The random potential is taken to be uncorrelated, normally 
distributed IISII with 



l^(x, r)y (y, T') = 2A 5(x - y)5(r - r') (2) 

where the overbar stands for the averaging over the disorder. The disorder is quenched so 
that the average of In W , where W is the appropriate partition function, is needed. 
The overlap for the above system in the continuum limit can be precisely defined as 

q(t) = \f^dT5{^,,{T)). (3) 

It can be obtained from the relation q{t) = —f '^'^^j^"'*'' \vo=o where /2(fo, t) is the free energy 
for the Hamiltonian in Eq.|l|. A scaling form 

f2{vo,t)=t^/^f{vot-^/n (4) 

is expected with determining the crossover exponent. This implies q = t^Q{vQt~'^^^) where 

^ = {^-<P-z)/z (5) 

This particular problem in a discrete version d = 1 has been studied by Mezard nu- 
merically in 0. His simulation results are consistent with C, = 1/z = 2/3, x = 1/2, and 
(f)/z = —2/3. A corollary is that the behavior of one chain remains unaffected by the 
presence of the other. This need not be surprising because the "screening" type effects in 
ordinary self and mutually avoiding polymers are finite density phenomena. Since at d = 1 
X — (j) — z = 0, one obtains ^(fo, t) ~ q{vot~'^^^). One of our aims is to determine 0. 
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Here we use the continuum formulation and map the problem to a KPZ type nonlinear 
differential equation for the free energy ^,17]. A dynamic renormalization group approach, 



a la Ref. [0, following a perturbative calculation in a Fourier conjugate space, is developed 
to study this equation. In the process, we obtain the scaling exponent for the interaction 
and establish that in any arbitrary dimensions x ^i-iid z remain the same as those of the single 
chain problem, as one would expect, even in the presence of the interaction. A series for the 
renormalized coupling constant can be identified by collecting the appropriate terms from 
the perturbative series. The recursion relation for the coupling constant, found after the use 
of the momentum shell technique, manifests the scaling form of the mutual interaction and 
the overlap. The exponents obtained through this process are for the critical point T^. for 
d > 2 but for the finite temperature phase for d = 1. An appeal to finite size scaling then 
enables us to extend the result to T ^ Tc for d > 2. 

We consider two chains which are tied at one end {t = 0) at the origin of the d dimensional 
space and extended upto xi and X2 at length t. The partition function iy(xi,X2,t), which 
is basically a sum of the Boltzmann weights of all configurations of two such chains can be 
written in the path integral form as 

iy(xi, X2, t) = Z'^"""''*^ 2)x; Px'2 exp [-n], (6) 

"'(0,0,0) 

where / Vx.2 stands for all possible paths of the two polymers and Ti. is given by Eq. 
|I]. This implies that the partition function satisfies a Schrodinger type equation written 
suppressing the argument of W 



7 + 17-^o(xi,X2,t) 

i=l,2 ^7 



W, (7) 



where (^olxi, X2, t) = V(xi, t) + V(x2, t) + fo(5(xi2) appear as the potential. 

Our approach starts with another version of the above equation for the free energy 
/;,(xi,X2,t) = (27/A)ln iy(xi,X2,t) which satisfies, again suppressing the arguments, 

|^=.E[7V?^ + ^(V./.n+^7o. (8) 

The impressive feature of this equation is that the parameter A, which was previously con- 
trolling the random potential and the mutual interaction between the chains in the original 
hamiltonian in Eq.|l], now appears only with the nonlinear term. Eq. ^ can be decoupled 
when there is no mutual interaction between the chains and such a decoupled equation can 
be solved exactly when A = |T^. One can then make a small perturbation in the non- 
linearity (i.e. in A which is effectively equivalent to introducing a small disorder into the 
problem. The recursion relation for vq, obtained after perturbation in the nonlinearity A, 
gives the influence of disorder on vq. 

A glance at Eq.^ shows that under the transformation x ^ b x and t —>■ t the 
parameters of the equation change as 

(7, A, ^o) ^ (6^-'7, b^'-^-'X, b^-'-^ vo). (9) 
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Therefore in the absence of the nonhnearity i.e when X = 0, z = zq = 2 and x = 2 — (i keep 
7 and Vq invariant. This scahng furthermore ensures the speciahty of d = dc = 2 since for 
d < dc a. small amount of nonlinearity becomes relevant with the growth of the length scale. 
Such speciality of dc is reflected later in the perturbative series. At this level one finds the 
necessity of going beyond the simple scaling analysis to carry out the RG analysis since the 
scaling dimension of fo, z — x — d, vanishes at ci = 1, while numerically it is found to be — 1 

§• 

The formal solution of Eq.H in (K,k, cu) space, Fourier conjugate to (xi,X2,t), is given 

by 

h{K,k,uj) = G'o(K,k,cj)^o(K,k,cu) - A/2 Go(K,k,cj) x 
Ip,,,n[P ■ (K - P) + q ■ (k - q)]h{p, q, Q)h{K -p,k-q,uj-Q), (10) 

where Go{'K,'k,u) = (7(K^ + k^) — iu)^^ represents the bare propagator and Jp^g^Q = 
J |^(^^- The fact that the random potential and the interaction are in the same foot- 
ing in the above equation is now utilised in defining the effective propagator G(K, k, u) and 
the effective coupling constant v as 

h{K, k, uj) = G{K, k, uj) [V{K, uj)6{k) + V{k, uj)6{K) + v6{K + k)6{uj)]. (11) 

Note that such a restriction on momenta is imposed automatically by the RHS of Eq.p!0|. 
Now we are in a position to initiate the perturbative series, the terms of which after 



disorder averaging leads to "closed loop diagrams". [See, e.g., Ref. fTj]. Here we shall con- 
sider terms upto O(A^) and 0{vo). Collecting the appropriate terms from the perturbative 
series for the renormalized propagator satisfying either K = or k = 0, one obtains a series 
identical to the renormalized propagator for a single chain in a random medium PJT7[]. It 
becomes evident from the series that there is no contribution at 0(A) since such terms either 
donot have any loop or even if they do, they contribute to 0{vq), a higher order term which 
we are not considering here. For convenience, we cite the series obtained for a single chain 
propagator with momentum variable K 

G(K, k, Lu) = Go(K, k, uj) + C(-A/2)2G2(K, k, u;){2A) q • (K - q) q • K x 

Go(K - q, k, - n)Go{q, 0, n)Go{-q, 0, -^]) (12) 

where C = 4 counts all the possible ways of noise contraction and k = in this case. The 
terms which contribute to the renormalization of the vertex follow the constraint K + k = 0. 
One can write down the series for vG(K, — K, 0) to obtain an equation in which the vertex 
and the propagator G(K, — K, 0) renormalization take place in a combined fashion and 
interestingly can also be isolated. To O(A^) the series for the renormalized propagator 
G(K, — K, 0) is same as Eq.^ with C = 8 and k = — K and v is given by the series 

v = Vo + 8(-|)2(2t;oA) ^^.(q ■ (K - q))^Goiq, -K, ^])Go(K - q, 0, -Q) X 

G'o(q,-q,0)Go(0,q-K,fi). (13) 

In a digrammatic representation the second term on the RHS of this equation would corre- 
spond to an exchange type diagram. At the face value, the series of G(K, — K, 0) does not 
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resemble the series for the single chain propagator 0(1^,0, u), but it is easy to show that 
there is no change in the renormalization of the line tension 7 from that of Ref. |^ in the 
long wavelength limit. This, furthermore, confirms that in any arbitrary dimensions, so far 
as the free energy, lateral extension are concerned, exponents do not change even after the 
inclusion of the mutual interaction. 

The second term of Eq.|l^ involves a momentum integration with an upper cutoff A which 
indicates the resolution upto which the system is probed. To study the variation of v with 
length scale, we execute an RG procedure consisting of two steps (i) an integration over the 
momentum shell between the momenta A and Aexp {—SI) and (ii) the momentum rescaling 
k —>■ kexp (—61) which restores the upper cutoff A as before. Now after carrying out the 
first step and the integration over Q we obtain an effective coupling constant (in the long 
wavelength limit) differing from vq by a term KdX'^voSl/2. Here Kd = with as the 

surface area of the unit d dimensional sphere and = is dimensionless. This additional 
term is the fluctuation contribution of the disorder and crucial for the RG analysis. A little 
manipulation after the rescaling (step (ii)) produces the recursion relation 

- = (z-X-c^). + ^, (14) 
with V = Vq when / = 0. The recursion for A is quoted below from Ref. 

which has a fixed point A* = [2d{2 - d)/{{3 - 2d)Kd)Y'^. At ci = 1 it is a stable fixed point 
which shows that the phase is influenced by the disorder at all temperatures. There is no 
physical fixed point for 1.5 < d <2. The nontrivial fixed point becomes unstable for d > 2. 
From the flow, one concludes that this unstable fixed point corresponds to a critical point 
that seperates the two phases dominated by the entropy (high temperature phase) or by 
the disorder (low temperature phase). The situation at c? = 2 is more complicated since 
disorder is marginally relevent [ll7| , p!9| . Using the relation x + z = 2, which is a consequence 
of Galelian invariance, and A = A* we find from Eq.[T^ 

^= {2 + d-2z) +d{d~2)/{3-2d). (16) 

At d = 1 exact values z = 3/2 and x = 1/2 yield the scaled variable as vt"^^^. The remarkable 
feature is that in Eq.|14| the term z — x — d which originates from the simple scaling analysis 
of Eq.|^, vanishes at d = 1 and the entire t dependence of the scaled variable comes only 
from the fluctuation part. The scaling exponent of v matches with the numerical prediction 
of Mezard at ci = 1 g]. Now, using the exponent 0, Eq.|ly we find 



J: = d{d-l)/[z{3-2d)] (17) 

Going back to the original problem of a single DP, the overlap is obtained by setting vq = 0. 
We therefore obtain 

q ~ t^, (18) 
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where, as mentioned before, this is the overlap in the low temperature phase for d < 1.5 but 
at the critical point for d > 2. The exponent S vanishes for d = 1, as has been found in 
Ref. 0. Since S < for d > 2, the overlap disappears at the transition temperature when 
the thermodynamic limit is approached - as one should expect. 

We extend our T = Tc result to the critical region by invoking a finite size scaling 
hypothesis, ||2^ q = C,^^Q{t/C,) where ^ is the longitudinal (parallel to t) correlation length 



and ^ ~| T — Tc I near the critcal point. Comparing with Eq.|TB], we find Si = S, so that, 
in the thermodynamic limit 

g-lT-Tj-'^^ (19) 

as T — Tc. Unfortunately the value of u is still not known with confidence ^j^- 

It is possible to explain the effect of random environment more physically. In the situation 
where the disorder dominates the physics, the chain is swollen to take advantage of the 
occasional traps that lower the energy. The loss in the entropy is offset by the gain in the 
energy, yielding ( > 1/2. In this scenario, it is therefore expected that the repulsion with 
another chain will have no significant effect. Hence to the leading order, vq is not to have 
any effect on the renormalization of the properties of the polymers. However, with two 
chains there will be a certain amount of overlap in their attempts to take the advantage of 
the same traps. On a bigger scale, such closely spaced traps would appear as an inetraction 
between the chains. This leads to the renormalization of vq changing its scaling exponent. 

In summary, we have shown through one loop RG analysis that the behavior of a sin- 
gle chain in the random medium remains unaffected even if we introduce another chain 
interacting with it through a short-range repulsive interaction. At d = 1 our results show 
a finite overlap at all temperatures indicating a strong coupling phase and the exponents 
match exactly with Mezard's numerical simulations. It may not be surprising that the 1 
loop result gives exact results in 1 dimension, because it is known to happen for all the 
exponents for the d = 1 DP problem. This is mainly a conseqence of Galelian invariance 
and fluctuation- dissipation theorem. In higher dimensions {d > 2), we evaluated the scaling 
exponent of the overlap at T = Tc. By using a finite size scaling ansatz, we thereafter predict 
the temperature dependence of the overlap as T Tc. 

I thank S. M. Bhattacharjee for fruitful discussions and comments on the manuscript. I 
also thank D. Dhar for several suggestions. 
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